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Abstract: The optimization of the coated metallic nanoparticles and 
nanoshells is a current challenge for biological applications, especially for 
cancer photothermal therapy, considering both the continuous improvement 
of their fabrication and the increasing requirement of efficiency. The 
efficiency of the coupling between illumination with such nanostructures 
for burning purposes depends unevenly on their geometrical parameters 
(radius, thickness of the shell) and material parameters (permittivities which 
depend on the illumination wavelength). Through a Monte-Carlo method, 
we propose a numerical study of such nanodevice, to evaluate tolerances (or 
uncertainty) on these parameters, given a threshold of efficiency, to facilitate 
the design of nanoparticles. The results could help to focus on the relevant 
parameters of the engineering process for which the absorbed energy is the 
most dependant. The Monte-Carlo method confirms that the best burning 
efficiency are obtained for hollow nanospheres and exhibit the sensitivity of 
the absorbed electromagnetic energy as a function of each parameter. The 
proposed method is general and could be applied in design and development 
of new embedded coated nanomaterials used in biomedicine applications. 
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1. Introduction 

Extensive theoretical and experimental studies of the interaction between an incident wave 
with metallic nanostructures, such as nanoparticles, nanodots [1], nanorings [2], nanocubes [3], 
nanoshells [4], have been achieved in physics and chemistry. Under illumination, these struc- 
tures, due to local plasmon resonances, are known to exhibit a high enhancement of the elec- 
tromagnetic field at their surfaces. This field strength depends on the size parameters of the 
particle and on two related quantities: the wavelength and the permittivity of the used materi- 
als (in the case of nanoshell, the core could be made of silica or just a vacuum). Due to this 
structural tunability of the plasmon resonances, such nanostructures can be of interest in a wide 
range of applications in biomedicine [5-7]. 

Therefore, these researches permit to develop applications in nanoscale plasmon waveg- 
uide [8] and nanosensing modalities [9-1 1]. In the past, the choice of the metallic materials and 
synthesized alloys, size and shape geometries of materials were studied experimentally [12-14] 
and interpreted as plasmon hybridization [15]. Therefore, the development of metallic nanode- 
vices, where the relative role of electromagnetic and resonant enhancement response could be 
precisely delineated, necessitates an accurate control of the local electromagnetic field enhance- 
ment near the metallic surface [16, 17]. 

For cancer photothermal therapy, nanodevices like coated metallic nanoparticles and 
nanoshells are currently used to burn cancer cells. Indeed, most biological tissues have a rel- 
atively low light absorption coefficient in the visible (VIS) and near infrared (NIR) regions 
(600-1300 nm) known as the tissue optical window or therapeutic window. Over this window, 
organic molecules have limited absorption [18], whereas gold nanodevices, which are biocom- 
patible nontoxic and easily conjugated to antibodies, absorb light millions of times stronger 
than the organic molecules. Then, almost all the absorbed light is converted to heat via series 
of nonradiative processes [18] to burn cells in which they are embedded. The plasmon reso- 
nance tuning helps to increase the absorbed light and therefore, the nanodevice capability to 
burn diseased cells. 

Until recently, in the absence of optimization of such devices, the absorption efficiency (i.e. 
the rate of absorbed energy by the particle relatively to the incoming illumination intensity and 
to the surface of the particle) is only of one order of magnitude [19,20]. In a recent study, 
it has been shown that it becomes possible to obtain an absorption efficiency of two order of 
magnitude in the most convenient wavelength domain for biomedical applications (around 900 
nm) [21]. Such a design optimization permits to obtain an increase in the absorption up to 
two order of magnitude in comparison with the classical coated nanoparticles. However, the 
sensitivity of such devices to manufacturing uncertainties has never been investigated. 

In this study, we focus on the sensitivity design of nanoshells used in photothermal therapy 
(i.e. the absorption of these particles). The main principle of the method is to select the particles 
parameters that keep the absorption efficiency greater than a threshold. These parameters are 
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within a range that defines the manufacturing tolerance. 

The paper is organized as follows. Section 2 is devoted to the description of the absorption 
efficiency computation and the requirements on absorption efficiency computation are defined. 
Then, in section 3, the optimisation scheme by the Monte-Carlo method is presented and ex- 
plained. The numerical setting related to the considered biomedical application (photothermal 
therapy for deep cancer using nanoshells) and results will be given and discussed in section 4. 
Finally, section 5 concludes the paper. 

2. Absorption efficiency computation 

The total electromagnetic field (near-field and far-field) properties of a metallic nanoobject de- 
pends intrinsically on the geometry and on the optical properties of the involved materials. 
Nanoshells are composed of a core with radius r\ and of a metallic coating or shell of thickness 
e (see Fig. 1 ). The core could be made of silica or just vacuum (the commonly used name of such 
nanoshells is hollow nanospheres with « i = 1 .0) whereas the shell is made of gold. The per- 
mittivity of the core and the permittivity of the coating are denoted £i and £2, respectively. The 
metallic material is characterized by a complex permittivity £2 = £ r + j£j, (where y 2 =-l). The 
nanoparticles are embedded in a non-absorbing medium with permittivity £„,, corresponding 
to biological surrounding. From these parameters (size and permittivities), the electromagnetic 
field can be computed accurately by a variety of methods, such as analytic Mie scattering the- 
ory for spherical geometries [22-26] and numerical methods such as Finite Element Methods 
(FEM), especially for nanoscale objects of more complex geometries [27-30]. In the following, 




Fig. 1. Nanoshell: inner radius r\ and shell thickness e. 

we compute the absorption efficiency Qafa (see Eq. (3)), which must be maximum to get the 
more efficient thermal effect. This absorption efficiency can be deduced from W a b s which is the 
rate at which energy is absorbed by the sphere [22]: 




where Ej, Hj, E s and H s are the incident and scattered electric and magnetic components, 
respectively. The integration is achieved on the solid angle dQ. = p 2 sin0d0d(j). It follows that 
the absorption cross section is defined by: 

C abs = ^ = ~%(2n+l){3i[a n +b n }- 

1 n= 1 



■\K\ 2 } 



(2) 
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where is the incident intensity, k = 2n^/e^,/X being the wave vector in the surrounding 
medium, and a„ and b„ are the scattering coefficients. The absorption efficiency Q u /, s is there- 
fore the absorption cross section C a y s per unit of area S = n(r\ +e) 2 : 

Qabs = %*= V27 2 , 2 Y J {2n + \){^[a n +b n ]-\\a n \ 2 + \b n \ 2 \\. (3) 

In the limit of small particles (/2/A << 1) [22], the computation of Q u /, s is often achieved by 
neglecting the order n > 1 in the expansion of the series (Eq. (3)) [22]. Therefore, in the limit 
of small particle approximation, the dipolar approximation consists in considering only the a\ 
and b\ (called dipolar electric and magnetic) terms of the series in Eq. (3). In the present study, 
the question is: is this dipolar approximation valid for the considered nanoshells (i.e. limiting 
the Q a i, s computation with n=l)? 

We show in Fig. 2 that even for small particles as nanoshells, this dipolar approximation is 
not valid. We summary in Table 1 the acceptable intervals of parameters for inner radius r\, 
thickness e, illuminating wavelength A and the relative real and imaginary parts of the bulk 
permittivity of gold 62(A) [31]. Figures 2(a) and 2(b) show the relative error between the ab- 



Table 1. Summary of Acceptable Intervals of Parameters r\, e, A, £ r (X) and £ f (A) 
Parameters 

inner radius r\ (nm) [1.0; 150.0] 
shell thickness e (nm) [1.0; 50.0] 

wavelength A (nm) [800; 1000] 

e r (A) [-42.0; -23.0] 
e,(A) [1.5; 3.0] 



sorption efficiency Q c ,b s computed from the dipolar approximation and the series (Eq. (3)), with 
60 terms, ensuring a convergence of the series, better than 10~ 12 ). The ratios of the radius to the 
wavelength are r 2 /X = 0.092 and r x /X = 0.089 (Fig. 2(a)), and r 2 /X = 0.022 and r x jX = 0.021 
(Fig. 2(b)). Even if the computation is based on the systematic variation of the permittivity of 




(a) (b) 



Fig. 2. Relative error between the exact computation of Q a b s and its approximation in 
the small particle limit for radii, as a function of the real e r and imaginary part £/ of the 
material permittivity for particle radii: (a) n /A = 0.089 and e/X = 0.003, (b) n /A = 0.021 
and e/X =0.001. 



the gold shell within the domain of wavelengthes of [800; 1000] nm, the external radius of the 
nanoshell is between 800 x 0.092 = 73.6 nm and 1000 x 0.092 = 92 nm and the relative er- 
ror is greater than 25% (see Fig. 2(a)). The relative error cannot be neglected in both cases. 
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This result shows that the computation of Q a j, s with dipolar approximation is not appropriate 
even for small particles with thin gold coating (see Fig. 2(b)). In the case of the present study, 
the series of the efficiency cannot be limited to the first order and the dipolar approximation 
is not relevant. The dipolar approximation cannot be used for the prediction of the location of 
plasmon resonance and therefore for optimization in the biological window [32]. 

Therefore, for sensitivity design, the multipolar computation of Q a i, s must be achieved (i.e. 
more than one term must be computed in the series). In the following, the computations of 
the <2afa are achieved by using the full Mie's theory (i.e. the series of Eq. (3) are computed 
with 60 terms, ensuring a convergence of the series, better than 10~ 12 ). In the next section, 
we introduce the computational Monte-Carlo method used to compute the sensitivity of the 
absorption efficiency Qab s to the various parameters of the problems. The tolerances for each 
parameter will be deduced, so that Q a j, s remains above a given threshold. 

3. The adaptive Monte-Carlo method 

The main purpose is to determine the parameters of the nanoparticles that are critical for med- 
ical application. For this, we introduce a boundary adaptive method based on a Monte-Carlo 
scheme. The proposed approach differs strongly from the metaheuristic optimization algorithms 
we applied in the previous studies, either to the inverse problem solving in near-field optics [33], 
or to optimize the Surface Plasmon Resonance planar biosensor [17, 34-37]. In those papers, 
the goal was to compute the best set of parameters to get the most efficient biosensor. Nev- 
ertheless, no propagation of uncertainty has been determined, despite the known uncertainties 
on permittivities and geometrical parameters. Recently, we have proposed a method to design 
the nanoshells with silica core for imaging and burning applications [21]. In that study, an 
evolutionary method (i.e. developping similar operational principles based on the evolution of 
the searched parameters as objective variables to reach a target) was used to optimize the size 
parameters and the permittivity of the shell using existing permittivities. However, the optical 
index of the core was fixed to that of silica (which is not the case in this study). The goal was to 
obtain at least 80% of the mathematical maximum of the absorption Q u \ os (resp. scattering) effi- 
ciency, for the burning (resp. imaging) purpose. This method is now used to find the best set of 
parameters to reach max{Q u bs) without fixing the optical index of the core nor the wavelength. 

Otherwise, in the present paper, the proposed method consists in fixing a tolerance on the 
optimum of the absorption efficiency, to determine a class of acceptable parameters. Then, the 
tolerance on the various parameters can be deduced and interpreted either in terms of sensitivity 
of the model or acceptable uncertainty in the process of fabrication. For this, a family of param- 
eters of the nanoshells (size and permittivities) is randomly generated, and the algorithm retains 
the maximal and minimal values (i.e. the boundaries) of the parameters for which the efficiency 
remains greater than a given fraction a of the optimal efficiency. The process is repeated (by 
regenerating randomly another set of nanoshells given the new boundary constraints and look- 
ing for additional "good" nanoshells) until a given number of "good" nanoshells are obtained 
and analyzed to depict the sensitivity to each size or optical parameter. This approach is derived 
from a classical metaheuristic engineering approach of the study of uncertainties propagation 
through a numerical model of the system, to compute the engineering tolerances. Actually, the 
Monte-Carlo method relies on repeated random sampling to compute their results (e.g. proba- 
bility of fitting some conditions) in a fixed domain (or boundaries) given usually by physical 
constraints. It is known to be an effective substitute to the traditional method of experimental 
plan, indeed to solve inverse problems [38]. 

In this paper, we present the method permitting to compute the intervals of optimized size 
parameters and core index that guarantee absorption efficiency Q a b s , given an acceptable per- 
centage a of its maximum. This maximum absorption efficiency max(Q a h s ) is corresponding 
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to the theoretical maximum best absorption efficiency obtained for the optimized size and ma- 
terials parameters. This maximum is computed with Eq. (3) and the best parameters are deter- 
mined through the evolutionnary method [21]. Therefore, the target of the problem is finding 
a family p of the optical and of the size parameters for the nanoshell, that enables to maintain 
Gnfes(p) > OCmax(Qab s ). The family of parameters is p = [r\,e, A, £1,62 (A)], where r\ is the ra- 
dius of the core, e is the thickness of the metal coating (shell), A is the illumination wavelength, 
£1 the permittivity of the core and £2 (A) the permittivity of the gold shell. The permittivity of 
the external medium is £ m . 

3.1. Principle of the method of model sensitivity study 

This numerical method is summarized as following, considering a threshold a = 90%, delin- 
eating an acceptable value of Q a b s : 

1. Define an initial domain of possible parameters (size parameters r® and e° and per- 
mittivities £j, ^(A)) by setting the boundaries according to physical constraints. 
The number of parameters set npt is fixed as the maximum of the domain of 
variation of the initial parameters divided by their uncertainty u in the process of 
fabrication [38]: npt = max{[max(rj) — min(rj ) )]/w(ri); [max(e°) — min(e°)]/w(e); 
[max(«j) — mm(n®)]/u(n\); [max(A 0 ) — min (A°)]/k(A)}. Within the range of param- 
eters, the maximum of Q a b s is computed with the above mentioned evolutionary 
method [21]. 

2. Random generation of a family of npt unknown parameters sets p(f ), using uniform laws 
in the considered domain, t is the iteration index initially set to 1 . 

3. Computation of the optimum geometrical parameters with Nelder-Mead method: this 
method requires an initial values for the optimum search. These initial values are 
the parameters set that gives the maximum value max ( (Qab s )> within the family: 
\r\ ,e, A, £2(A)], £1 being unchanged. This maximum is voaxj^ti{Qabs)t an d varies slightly 
with t . Actually, the basic Nelder-Mead method is a multidimensional unconstrained non- 
linear method that cannot handle with boundaries in the search domain and therefore, if 
n\ would be also a parameter to be optimized, the output of the algorithm would give 
n\ =0, which is not compatible with the physical constraint on the optical index of a 
dielectric medium n\ > 1. The optimum set generated by this method will be denoted 
Pmax(t) and the corresponding absorption efficiency maxNM(Qabs)- 

4. Selection and storage of the best parameters sets of the family p(f) for which Q c ,b s > 
a maxNM(Qabs)- The retained sets form a family denoted p a (f)- The set p max (t) is not 
added to this family p a (t) but is used only to determine max NM (Q a b s ). 

5. Updating boundaries of the domain using the minimum and the maximum of each pa- 
rameter in p«(f )■ 

6. Increment t and loop on step 2, until the size of all the selected sets of parameters |J; Pa (0 
is greater than npt . The final number of iteration is N t . 

The associated tolerance in fabrication will be deduced from the boundaries of this last family. 

3.2. Benchmark of the adaptive Monte-Carlo sensitivity study 

At this stage, two points should be clarified. First, the advantage of the adaptive method v.s. 
the classical Monte-Carlo method. To reach an accuracy on the boundaries determination lower 
than 1%, the number of random parameters should be greater than 100,000, the convergence of 
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the classical Monte-Carlo method being in 1 j y/npt. The advantage of our adaptive method lies 
in the adaptation of the domain of search, keeping constant the number of random parameters 
at each iteration, from a generation to the following. Therefore the convergence is better at each 
new iteration and the precision on the boundaries is increased. 

Second, the numerical cost of the method has to be compared to the cost of the four 
loops required in a systematic study. For this, we use the numerical parameters, that are 
physically acceptable for the initialization of the algorithm. We consider r\ E [1 ; 100] = 
[min(rj);max(rj)] nm; e £ [1;50] nm as small particles are preferred in burning applications. 

For the benchmark of the algorithm, we used the same parameters as in Ref. [21] where 
non heuristic method was used to compute the tolerance on the parameters. Table 2 shows the 
various parameters of the study and the obtained results. In that study the optical index of the 
silica core (resp. external medium) was £i = 2.4 (resp. £3 = 1.77). For illustration, we focus 
on a specific result for gold nanoshell and A = 850 nm. The threshold was a = 80%. The 

Table 2. Benchmark of the Adaptive Monte-Carlo Model Sensitivity Study: Comparison 
with Systematic Study [21] 



Parameters 


Ref. [21] 


Monte-Carlo 


inner radius r\ (nm) 


25.1 ±9.8 


[16.7(0.6);30.2(0.6) 


shell thickness e (nm) 


3.0± 1.3 


[1.9(0.1);3.9(0.1)] 


max NM (Q abs ) 


11.6 


11.6 


best r\ (nm) 


25.1 


24.7 


best e (nm) 


3.0 


2.9 


Number of Q a f, s evaluations 


24,000 


6700 



*The maximum of the efficiency is deduced from Table 2 in Ref. [21] and the reference by Loo et al. [5]: 
1 1 .6 = 0. 144 x 80.2. The number of evaluations is that required for permittivity choice as well as optimization of 
the geometry of the nanoshell: 20,000 + 4,000 [21]. For the Monte-Carlo method, the standard deviation of the 
boundaries across the iterations is also indicated between parenthesis. 

computation of the standard deviation of the boundaries obtained for all iterations can be held 
as an indicator of the confidence on the computed boundaries. Indeed, the standard deviation of 
the boundaries can be considered as uncertainty and is an indicator for limiting the significative 
digits of the results. Increasing npt results in a 1/ 'y/npt decreasing of the standard deviations 
of the boundaries, the number of iterations being almost constant (N t around 40). 

The Monte-Carlo method, gives the tolerance on the radius r\ and the thickness e. These 
tolerances are deduced from the intervals of acceptable values of each parameter. The tolerances 
correspond to those obtained from a systematic study like experience plane [21], but the number 
of required evaluations of Q a b s is reduced. The computational scheme requires less than 40 
iterations (with 990 evaluations of Q a ts) an d therefore, less than 40,000 evaluations of Q a b s , 
while a three loops systematic study would require 990 x 49 x 200 ~ x 10 7 evaluations [21]. 

These results are validating the Monte-Carlo approach which helps to compute sensitivity of 
the model to the input parameters. From Table 2, we can deduce that the critical parameter is the 
thickness of gold, the tolerance on r\ being around four times greater. After convergence of the 
Monte-Carlo scheme, the boundaries of the last hypercube can be used to define the tolerance 
on each parameter. 

A last bench has been made before the use of this method for more general problems. Indeed, 
it is well known that the reiteration of metaheuristic methods has to be made to characterize sta- 
bility and dispersion of the results. Thousand realizations of the proposed Monte -Carlo code, 
with different initializations of the random generator, have been made to check the stability of 
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the algorithm. The obtained results show a standard deviation of the boundaries, over the thou- 
sand realizations, which is smaller than the initial tolerance on each parameter. Therefore, the 
method could be applied with sufficient confidence to the sensitivity design of nanoshells. Let 
us investigate the case of hollow nanoshells {n\ 1) and coated silica (n\ = 1.54) nanospheres 
in fat instead of water [39]. 

4. Sensitivity design of nanoshells for photothermal therapy 

4.1. Sensitivity of the model with free optical index of the core 

We focus on the design of gold coated nanoparticles (or nanoshells) for applications in the 
biomedical domain of wavelengths (A G [800; 1000] nm), for cell burning purpose. The control 
of absorption efficiency by tuning the wavelength and the size in visible and infrared region has 
been proved experimentally [19, 34]. Nevertheless, the optical constant of gold nanoparticles 
is not known exactly. Moreover, strong variations have been found for bulk and obviously for 
thin layers, depending on the mode of deposition. With such variations in the values of the 
permittivity for gold, the question is: is the numerical optimization using the bulk permittivity 
for the gold shell adapted? In fact, despite the variations between values of the permittivity for 
gold nanoparticles and bulk values, the results of computation of cross sections are in good 
agreements with experiment results, especially in the case of spectroscopic studies of such 
nano-devices [40-42]. Hopefully the following results confirm that the geometrical parameters 
are more critical in the design than the optical index of gold. 

The permittivity of the core E\ G [1 ; 16] as the commonly used core is made of silica or 
vacuum (with vacuum core, the corresponding nanoshells are commonly known as hollow 
nanospheres). The uncertainty on the optical core index n\ = ^/e[ is fixed to 0.01. The permit- 
tivity of the gold coating (or shell) is obtained from the Palik's data [31], in the [800; 1000] nm 
range of wavelength A: £2 (A). Therefore, the real part of the gold permittivity lays between 
—42 and —24. The imaginary part of this permittivity is between 1.5 and 3. A fit of this permit- 
tivity with A's step equal to 1 nm is used. On the contrary of Ref. [21], the material of the shell 
is fixed and our results reveal the shallow sensitivity on this parameter. 

Let us investigate first the most general case, where the index of the core must 
also be determined, the initial set of parameters is summarized in Table 3. Therefore, 

Table 3. Parametric Setting: Domain and Accuracy 

Parameters Domain uncertainty u 

inner radius r\ (nm) [1; 100] 0.1 

shell thickness e (nm) [1;50] 0.1 

illumination wavelength A (nm) [800; 1000] 1 

core optical index « 1 = [L4] 0.01 

the size of each family of parameters at each iteration, npt, can be evaluated: npt = 
max{(100-l)/0.1,(50-l)/0.1,(1000-800)/l,(4-l)/0.01,(800-400)/l} =990. 

In the following, we present the histograms for each parameter (see Fig. 3). These results are 
obtained for a value of Q a t, s which falls within the interval [Q.9maxNM{Qubs)\ max NM(Qabs)]- 
At this stage, the determination of the tolerance on experimental parameters can be evaluated: 
it is actually the size of the last domain [38] (i.e. the minimum and maximum on the x-axis of 
each histogram). The boundaries of this domain can be considered as an absolute confidence 
interval of the parameters for the threshold a (all parameters within this interval product ac- 
ceptable efficiency). Nevertheless, the visual inspection of each histogram in Fig. 3 enables the 
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12.8 13 13.2 13.4 13.6 13.8 14 



(e) 

Fig. 3. Histograms of (a) the wavelength, (b) the optical index of the core {n\ = *j£\) (c) the 
radius of the core, (d) the thickness of the shell, (e) the absorption efficiency. The relative 
frequency is plotted in percents, and the number of class is deduced from the uncertainty 
in Table 3, except for the absorption efficiency where the size of each class is fixed to 
0.5%max NM (Q abs ). 
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classification of each parameter in terms of sensitivity of the result. The conclusions are: 

• hopefully, the sensitivity of Q a i, s to the wavelength (and consequently to the permittivity 
of gold) is low. In the investigated range of wavelengthes ([800; 1000] nm), the permittiv- 
ity has monotonic variations. Therefore the nanoshell can be efficient in this whole do- 
main of wavelengthes. Moreover, the exact knowledge of permittivity of the shell seems 
not to be critical. This behavior can be explained by the broadening of the multipolar 
resonance of the nanoshell due to the non negligible imaginary part of the permittivity. 

• The less the inner permittivity is, the better is the efficiency of the nanoshell (see Fig. 3(b) 
with n\ w 1.0 and Fig. 3(e) with 12.7 < Q a bs < 14.0). This confirms that the hollow 
nanosphere seems to be the more efficient. Schwartzberg et al. have shown that the ab- 
sorption band of these particles can be tuned by adjusting the thickness of the gold shell 
and the inner radius and thus would enable both strong scattering and absorption effi- 
ciency [34]. 

• The computed tolerances on radii are within these of fabrication processes (0.6 nm [34, 
35]). But the thickness of the gold shell is the most critical parameter: e = 2.1 ±0.6. r\ = 
22.5 ±6.5 and the best parameters are r\ = 22.2 nm, e = 2.0 nm, to get max^MiQabs) 

1 4.0. The best parameters for A, n and e are close the center of the intervals of tolerance 
and n\ w 1.0. The standard deviation of the boundaries of intervals are respectively 0.5 
for r\, 0.05 for e and 10~ 3 for n\. These values are around ten times lower than the size 
of each interval, ensure the validity of the above intervals of tolerance, and confirm the 
sensitivity of Q a j, s to each parameter. In particular, the high sensitivity of Q a j, s to the 
thickness of gold is observed and can be assimilated to the same behavior observed for 
Surface Plasmon Resonance based planar biosensors [17]. 

• The histogram of Q a b s shows a maximum of 14.0 and of course a minimum of 12.9 which 
is related to the threshold a = 0.9. 

With such parameters, the absorption efficiency, relatively to the reference particle size (see 
Loo et al. [19]), corresponds to a gain of 60 to 100. 

Figures 4(a)-4(d) show the convergence of the boundaries of the domain. In this case, the 
maximum number of iterations required to get at least npt parameters sets UrPa( f ) satisfying 
the target Q a b s {p G [0.9maxAfM(2af>.v); max A'M(2«fo)]) is N t = 30. It illustrates the rapid conver- 
gence of the boundaries which can be considered as the tolerance on each parameter, and their 
low dispersion. The minimum and the maximum of the wavelength (not shown here) remain 
the same along the iterations. Figures 4(a)^-(c) show the quick convergence of ri\ toward unity, 
and of the shell thickness in a small interval compared to the initial one. Clearly, the sensitivity 
of Q a bs to r\ is less. Figure 4(d) shows that the maximum of absorption among the family p a (f ) 
(diamonds) and the best solution p„ mx {t) determined by Nelder-Mead optimization (plus), are 
close together (max(Q a b s ) ~ 14) since iteration 20. This confirms that the space of parameters 
search is therefore reduced efficiently. 

4.2. Sensitivity of the model with silica core 

A second application is carried to study one of the commonly manufactured nanoparticles: the 
nanoshells with silica core. For this, the optical index n\ is fixed to 1.54. This is possible as 
its variation with the wavelength within the range [800; 1000] nm is negligible. The proposed 
Monte-Method results are r\ G [18.2(0.6);28.2(0.6)] nm and e G [2. 1(0.06);3. 3(0.06)] nm. 
The standard deviation of each boundary in the family Ur Pa(0 of 990 acceptable parameters 
is indicated in parenthesis. The tolerance on the geometrical parameter is of the same order 
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Fig. 4. Example of the convergence of the boundaries as a function of the iterations: (a) the 
optical index of the core (n\ = -y/Ei) (b) the radius of the core, (c) the thickness of the shell, 
(d) the absorption efficiency. The minimum and the maximum for each parameter for the 
last iteration are the final results. The plot of the convergence of the wavelength is useless 
since its interval remains almost the same at each iteration. 



of that in the previous subsection, and the best values are comparable. Therefore, the core 
index n\ of the nano-shell can be considered as much less critical than the thickness of gold e. 
These results are confirmed if we decrease the threshold of tolerance to a = 50%. In this case, 
n G [9(1);46(1)] nmandee [1.0(0.2);5.9(0.2)] nm. 

Moreover, max^MiQabs) = 12.2 for r\ = 23.2 nm and e = 2.7 nm and therefore, the best 
efficiency of the coated particle is hardly smaller than 0.9maxMv/(Gafo.s) for hollow spheres. 
This result shows that the gold coated spheres may continue to be used to burn cancer cells. 
This low loss of efficiency should not prevent their use for their lesser fragility. Nevertheless, 
the control of the thickness of the coating and its quality should be improved mainly as forming 
a uniform shell on the silica core is very difficult for small clusters [43]. 

5. Conclusion 

We proposed an iterative Monte-Carlo method with boundary adaptation, to compute the sen- 
sitivity of the absorption efficiency of nanoparticles. The goal was to determine the sensitivity 
of the geometrical and material parameters, that could maintain a sufficiently high absorption 
efficiency, to get an elevation of its temperature for burning purposes. The proposed algorithm 



#1421 19 - $15.00 USD Received 2 Feb 201 1; revised 22 Apr 201 1; accepted 6 May 201 1; published 17 May 201 1 
(C) 201 1 OSA 1 June 201 1 / Vol. 2, No. 6 / BIOMEDICAL OPTICS EXPRESS 1595 



enables to spare evaluations of the target (Qabs) with regards to a systematic study, through 
simple loops. The obtained results are coherent with former experimental and theoretical re- 
sults and are statistically significant. In particular, the method shows that the hollow nanoshells 
(no material in the core) are more efficient. The extent of the domain of possible parameters 
shows that the shell should be synthesized with more precision to guarantee a maximal absorp- 
tion. On the other hand, the efficiency of the nanoparticles also depends on the gold process 
of deposition and on the inner radius. The method introduced in this study is general and can 
be applied to a wide range of problems and especially to nanoparticles in order to obtain an 
optimized design of the size, the shape, the geometry and to select the most adapted for a spe- 
cific application. The extension to non-spherical particles with complex geometry may give 
perspectives for both research and applicative purposes, especially in the optical imaging and 
in the biomedical domain. The extension to complex or non-spherical geometries will neces- 
sitate the definition of additional geometrical parameters (like roughness), the development of 
an adapted numerical method (e.g. FEM, DDA,..) to compute the target to optimize (Qabs, or 
absorbed enery) for the considered geometry. Nevertheless, this study opens a way in the eas- 
ier manufacturing and development of a new family of sophisticate nanostructures (nanoshells, 
nanorings, nanodots, etc.) with optimized design. 
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